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Abstract. The scalar wave equation is considered in the background of a charged Vaidya 
metric in double null coordinates (u, v) describing a non-stationary charged black hole with 
varying mass m(v) and charge q(v). The resulting time-dependent quasinormal modes are 
presented and analyzed. We show, in particular, that it is possible to identify some signatures 
in the quasinormal frequencies from the creation of a naked singularity. 



1. Introduction 

Even though macroscopically charged black holes are not expected to be produced in realistic 
astrophysical scenarios, we can obtain fascinating information on the structure and behavior of 
Einstein equations and general relativity with the study of these objects. Here, we used the 
charged Vaidya metric to model time-dependent black holes with an electric charge, and study 
the behavior of the quasinormal modes in a non-stationary background. 

The charged Vaidya metric [1] [2j [3] is an asymptotic flat and spherically symmetric solution 
of Einstein-Maxwell equations in the presence of a radial flow of a null-type fluid. For our 
purposes, it is more convenient to consider the charged Vaidya metric in double- null coordinates, 
as proposed by Waugh and Lake in [4j. In double-null coordinates (u, v, 6, 4>), the charged Vaidya 
metric reads 

ds 2 = -2f(u,v)dudv + r 2 (u,v)dn 2 , (1) 

where dti? stands for the metric on the unity sphere (6,4>). The Einstein-Maxwell equations 
imply that the functions f(u,v) and r(u,v) and h(u,v) obey 

f(u,v) = 2B(v)r >u (u,v) (2) 

and 

/ \ / 2m(v) q 2 (v)\ , . 

r, v (u,v) = -B(v) (l - — U + , (3) 

where B(v), m{v) and q(v) are arbitrary integration functions. With the choice of B(v) = 
constant, one can interpret m(v) and q(v), respectively, as the total mass and the total charge 
of the solution. Both functions are only required to be monotonic, but a further condition must 
be imposed due to the weak energy condition in the case where both m and q are non null [5]. 

The solutions of Eq. ([3]) corresponds to the congruence of ingoing and outgoing null geodesies 
of the metric ([1]). We adopt here the semi-analytical method proposed in [6l [7] to construct 
such geodesies. In Fig. [T]we present as an example the null geodesies obtained for a case with 
m = q = 0, that is, a Reissner-Nordstrom (RN) black hole, in this parametrization. Fig. [2] 
presents the same solutions for a time dependent background. For an interesting analysis of the 
time dependent horizons in the Vaidya metric, see |8j. 
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Figure 1. Example of null- geodesies obtained 
from eq. ([3]) for a RN black hole, with B = 
-1/2, m = 0.5, q = 0.4, n = 4 and A = and 
taking r(u,v = 0) = —u/2. Here r + and r_ 
are defined as usual as r± = m ± i/m 2 — g 2 . 



Figure 2. Same as in Figure [TJ but this 
time for a time dependent charged black hole, 
with the choices B = —1/2, m(v) = 1.0 and 
q(v) = l(q 2 -qi)[l + tanhp(v-v m )]+q 1 , with 
qi = 0.5, c/2 = 0.7, p = 0.1 and v m = 90. 



2. Scalar perturbations 

The evolution of a massless scalar field in this background is given by the Klein-Gordon equation: 



1 



-9 



Decomposing the field ^f(u,v,9,(j)) in spherical harmonics 

£,m 



and using eqs. ([2]) and ([3]), we find 



with the potential V(u, v) given by 



f,uv + V(u,v)cp = 0. 



f ) = 2/^' U ) 1 r 2 + ^3 



(4) 



(5) 



(6) 



(7) 



For more details on the numerical algorithm, initial conditions and wave extraction process for 
Eq. ([6]) we refer the reader to refs. |£l HOj . 



3. Numerical results 

We first tested our code for static Schwarzschild and Reissner-Nordstrom configurations, with 
very good agreement with values found in the literature |11[ [T2] . 

In Fig. ([3]) we present the time variation in the uj r x uji plane of the real and imaginary parts 
of the QNM of the scalar perturbation with i = 2 for a RN black hole with constant mass and 
a time dependent charge function q(v). The frequencies change in a similar way with the time 
variation of the mass or the charge, but are typically more sensitive to the charge. The specific 
case depicted in Fig. ([3]) corresponds to a charge variation where q/m goes from 0.7 to 0.8, where 
we see how the variation goes from one stationary state to another through a non-adiabatic 
trajectory. The relevant timescale here for characterizing mass or charge changes which are 
adiabatic or not is essentially the same one of the q = case studied in |10] : QNM modes will 
relax and follow an adiabatic trajectory in the oj r x plane whenever r'!(v) < uJi(v) . 




Figure 3. Left plot: quasinormal modes in the u r x oji plane, showing the transition between 
stationary states for the i = 2 scalar perturbation of a black hole with constant mass m = 0.5 
and q(v) = \{q2 — + tanhp(w — v m )\ + ql, with q\ = 0.7, q2 = 0.8, v m = 75 and different 
values of p. Right plot: the minima of 0Ji obtained in the transition shown in the left plot for 
different values of p, to quantify the non-adiabatic behavior. 



We would like also to use our formalism and numerical setup to probe the formation of a 
naked singularity in the spacetime. In order to do so and avoid numerical problems as q — > m, 
we choose m = mo and q(v) = (q2 — <?i)[l + tanh/)(u — v m )] + qi, taking qi = tuq. This implies 
that the RN black hole will be (rigorously) extremal only for v — > oo, but allows us to probe 
configurations very near this end state at finite times. By varying the parameter p, we can 
control the speed with which the black approaches the extremal state. The results are shown in 
Fig. In the left plot we show how the results depend on the manner in which q(v) —> m. 
We remind here that an extremal RN black hole is not a naked singularity, and its QNMs are 
well defined as the appropriated limit of the non-extremal case |13} [1^1 [To"] . The results in the 
right plot of Fig. [¥] show what happens when we try to artificially create a naked singularity. We 
chose configurations with constant charge q and time dependent mass functions with final mass 
m,2 approaching and smaller than q. The results for oj r are shown in the right plot of Fig. U]and 
are truncated when m(v) = q since the simulation crashes at this point, as expected. 
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Figure 4. Left plot: same as in the left plot shown in Fig. (|3J), but this time with q^ = 0.5 so 
that the RN black hole becomes extreme as v — > oo. Right plot: co r (y) of a scalar perturbation 
with I = 2 obtained for "artificial" background configurations with m2 approaching and smaller 
than q; 0Ji(v) presents a similar behavior. 



4. Final remarks 

Our results confirm, for the charged case, the same non-stationary behavior corresponding to 
some sort of inertia for QNM frequencies identified in [9]. In particular, for situations with 
f'i-{v) > ^i(v), the QNM frequencies will not follow in the w r xw,; plane a trajectory corresponding 
to the instantaneous frequency associated to a RN black hole of a given q/m ratio, and a certain 
inertial behavior for the QNM can be identified. Moreover, we see that the faster the change in 
m(v) or q(v) is, the bigger the deviation from the stationary regime will be (see Fig. (3)). These 
points are now under investigation, and a more complete analysis will appear soon |16j . 

As for the q = case considered in [S] , an interesting extension of this work would be analysis 
of the highly damped QNM, the so-called overtones. Since for such overtones the ratio \ui/uj r \ 
is always larger than for the fundamental (n = 0) QNM considered here, including, for sufficient 
large n, cases for which \uJi/u) r \ > 1, it would be interesting to check if the non-stationary 
inertial behavior could somehow be attenuated for n > 0. We notice that the numerical analysis 
presented here cannot be extended directly to the n > case since one cannot identify the n > 
frequencies with sufficient accuracy. We believe this could be attained, in principle, by means of 
the WKB approximation. 
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